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Abstract 

We have performed self-consistent calculations for first and second row atoms 
using a variant of density-functional theory, the optimized effective potential 
method, with an approximation due to Krieger, Li and lafrate and a correlation- 
energy functional developed by Colle and Salvetti. The mean absolute deviation 
of first-row atomic ground-state energies from the exact non-relativistic values 
is 4.7 mH in our scheme, as compared to 4.5 mH in a recent configuration- 
interaction calculation. The proposed scheme is significantly more accurate than 
the conventional Kohn-Sham method while the numerical effort involved is about 
the same as for an ordinary Hartree-Fock calculation. 



1 Introduction 



Since the development of modern density-functional theory (DFT) by Hohenberg, 
Kohn and Sham |^ the accurate treatment of exchange and correlation (xc) 
poses a major challenge. Although there exist quite a number of approximations 
for the xc-energy functional [|, ||, ||, 0, ||, 3, which give rather good results 
[p!0| , pmi , the need for further improvement is undisputable. 

In conventional DFT the xc energy E^c is approximated by explicit functionals 
of the density p. The local density approximation (LDA) and the generalized 
gradient approximations (GGA) ||, |^ fall in this category. However, there 
exists another approach within the framework of DFT, the so-called optimized 
effective potential (OEP) method first suggested by Sharp and Horton [|1^] and 
refined by Talman and Shadwick [|13| . In this approach, the total energy (atomic 
units are used throughout) 

2 J J \r-r'\ 
+E.c[{Vj^}] (1) 

is a functional of = A^^ -|- A^j^ spin orbitals (pj^j resulting from a single-particle 
Schrodinger equation with a local effective potential: 




= ej„ipj^{r) j = 1,...,N^ cr =t, i . (2) 

The optimized effective potential, V^^^{r), is determined by the condition that 
its orbitals be the ones that minimize the energy functional (|1|). The stationarity 
condition 



OEP 



= (3) 

V=VOEP 



may be written by virtue of the chain rule for functional derivatives as 



E/^r' 



+ c.c. = 0. (4) 

y^yOEP 



Given an approximation of E'ajc [{v^jo-}], the first derivative in the integrand of 
equation (§) is easily computed from equation (P. To obtain a more accessible 
form of the second derivative in equation (H), one may use perturbation theory 
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for an infinitesimal perturbing potential 5V^(r) to get 



SVJr) 



E 

fc=i 



^ka 



(5) 



Insertion of the functional derivative of (|I|) and equation into equation 
leads, via the one-particle equation (Q), to the following integral equation: 



E / dr (K^r (r') 



( 
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\k^i 
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Via[r)^,A'^ ) + c.c. 



with 



and 
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We emphasize that from a fundamental point of view the total-energy func- 
tional (|lD is identical with the Hohenberg-Kohn energy functional 

EZ^[P^.Pi]=Ts[p^.Pl]-Z j ^-^dv + \j jP^0^dvdv' + E^^[p^,p,] (9) 

of conventional spin-DFT: By virtue of the Hohenberg-Kohn theorem for nonin- 
teracting systems, the spin orbitals ipja are functional of the spin densities po- so 
that 

E^onP^.Pi] = E?r[W,AP^.Pi]}]- (10) 

In other words, the orbital functional (|l|) is an implicit density functional. If 
the (unknown) exact xc functional were used in equation (|I]) then the optimized 
effective potential determined by equation would be the exact Kohn-Sham 
potential. Any approximation of E^c in equation (|ip, on the other hand, leads to 
an approximate Kohn-Sham potential. 

There are three non-trivial density functionals contributing to the Hohenberg- 
Kohn-total-energy functional @: the non-interacting kinetic energy functional 
T5[p|,pJ, the exchange part Ex[p],pi\ and the correlation part Ec[pi;,p\\ of 
Exc [PT'Pi]- 1^ "Es [P]iP\\ is approximated by an LDA, one obtains the Thomas- 
Fermi model. The transition from Thomas-Fermi to modern Kohn-Sham theory 
is equivalent to replacing the approximate LDA functional by the exact orbital 
representation 



t: 



exakt 



[PhPi] 



E E 



(pi„{v)dY. 



(11) 
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The transition from standard Kohn-Sham theory to the OEP method can be 
viewed in much the same way: While in ordinary Kohn-Sham theory the exchange- 
energy functional is approximated by LDA or GGA-type functionals such as the 
one due to Becke [§, the OEP method allows one to employ the exact orbital 
representation 

of the Kohn-Sham-exchange-energy functional^, which leads, by virtue of equa- 
tion d^), to the variationally best local exchange potential V^^^{r). Of particular 
importance is the fact that V^^^{r) is manifestly self-interaction free. Of course, 
the correlation part Ec [pti Pi] still has to be approximated, but even for Ec [p^, pj.] 
the representation in terms of orbitals allows more flexibility in the construction 
of approximate functionals. In particular, V^^^{r) is currently the only ap- 
proximate xc potential featuring the required discontinuities as a function of the 
particle number Ncr fl^ . 



The full OEP method, however, has a serious drawback: The solution of 
the integral equation (P) is numerically very involved. Recently, Krieger, Li 
and lafrate [T3, W\ have proposed an approximation (which we will refer to as 



KLI) to the OEP integral equation (^, in which V^^^{r) is essentially obtained 
from the solution of two linear [No- x N^j) equations. They also showed that 
in the x-only case, the results of the KLI approximation ||18|, |19|, ^ are nearly 



identical with those of the exact x-only OEP-method ||21| , P2| , |23|] . In addition, 
the KLI approximation preserves all of the important advantages of the exact 
OEP method, such as the correct asymptotic — ^ decay of V^^^{r) and the 
discontinuities of V^^^{r) as a function of the particle number. 

So far, the KLI approximation has been used exclusively in the x-only limit. 
In this letter, we present the first KLI calculations which include correlation 
effects. From the various available correlation-energy functionals we choose the 



one developed by CoUe and Salvetti [E^ pH|. This orbital-dependent functional 



gives correlation energies within a few percent of exact values if Hartree-Fock 
orbitals are inserted. The transformed version of it, the density functional of 

^The so-called Kohn-Sham-exchange-energy functional |23) is distinguished from the 
Hartree-Fock exchange energy expression by the fact that the orbitals in come from a 
local single particle potential. It is this restriction to local potentials that allows one to estab- 
lish the Hohenberg-Kohn 1-1 correspondence between densities and potentials which implies 
that the orbitals in equation ( p^ and thus and Ec are implicit functionals of the densities. 
If the restriction to local potentials is dropped, then only a 1-1 mapping between (non-local) 
potentials w{r,r') and one-particle-density matrices p(r,r') can be established |Q. In this 
formalism Hartree-Fock exchange is treated exactly but the correlation energy becomes a func- 
tional of the density matrix rather than the density, which can lead to a 7V-representability 
problem [^. For the density, on the other hand, A^-representability is always satisfied [p^ [T7[|. 
Therefore the restriction to local potentials and densities (rather than non-local potentials and 
density matrices) is essential. 
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Lee, Yang and Parr 0, has been shown to perform very well in self-consistent 
calculations p6|, especially in connection with an exchange energy functional 
proposed by Becke 



2 Theory 

In order to find an approximate solution to the OEP integral equation 



Krieger, Li and lafrate started from the Slater exchange potential given 
by 

Vi{v) = ^p,,(r)n„<,(r), (13) 

PoyT) i=i 

with 

P„(r) := |»Ji<,(r)P (14) 

«...(r) 1 ^^Mll, ,15) 

where E^a for cr =t, | represents the spin- up and spin-down contribution, respec- 
tively, to the right-hand side of equation (p!2D. KLI [1^ argue that a generalized 



form of this potential, which differs from equation (|r^) only through orbital de- 
pendent constants Ci„ added to Uxiaiji^), should yield a better approximation than 
Slater's formula. Moreover, they were able to show through a series of manip- 
ulations [0, that the exact V^^^^ijc) may be - up to within a small term they 



neglected - cast into a form similar to (pTSf ), namely 



1 

K'^r (r) ^ V;t^'(r) = ^ E h«.(r) + (l^J^^ " ^.c..)] (16) 

where the constants {VxcFJ' ~ ^xda) are the solution of the linear equation 

E (5,. - M,,,) iy^^f - u,,,^) = V,%, - j = l,...,K-l (17) 



4 = 1 

with 

M,..:=/^r^^^^%M, (18) 

TV / \ 

Kt(r):=E^«-^W- (19) 

Uxcja denotes the average value of 'U2;cicr(r) taken over the density of the ja orbital, 
i.e. 

Uxcja = J pjaiT^)Uxcja{r)dv (20) 
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and similarly for V^^^. The indices in equation ( ]T7| ) run over all occupied orbitals 
except the one corresponding to the highest single-particle energy eigenvalue 6^^^. 
If the highest occupied orbital is degenerate all orbitals with the energy e^a- are 
excluded from the linear equation ([T^) . 

Splitting the exchange-correlation energy E^c, as above, into an exchange and 
a correlation part, the exchange part (|1^) leads via equation (|^ to the following 
expression for u^ja- 



1 



r — r 



(21) 



The correlation part, E'c, within the approximation by Colle and Salvetti 



25 



E. 



is given by 
= - a6 y" 7(r)^(r) 



ab 



Vp(r) 



^Ep.(r)Ap.(r) + ^p(r)Ap(r) 



dr 



where 



-a / 7(r)44 
7]{r) 



7(r) 
r]{r) 



4 PT(r)Pi(r) 
p(r)2 ^ 

1 + dp{r)~3 
P(r) 



'3 6' 



-cp(r) 3 



r/(r) 



(22) 

(23) 
(24) 
(25) 



The constants a, b, c and d are given by 

a = 0.04918, 6 = 0.132, 
c = 0.2533, d = 0.349. 

Performing the functional derivative with respect to the one-particle orbitals, one 
obtains for ?/cjo-(i') 



ri{r 
ab 1 



7(r) + p(r) 



1 



d'j{r) 



4 vUr) 



d 



7(r)^(r; 



; 3 ^ r/(r)2 
4EP^'(r)ElV^-'W 



■(Vp(r))' + (pt(r)Api(r) +Pi(r)Apt(r) 



6 



ab 



- v(7(r)e(r))(Vp(r) + Vp,(r; 



A 7(r)e(r) /9^(r 



2 
ab 

T 



- ab 7(r)^(r) 



+ ab 



+ ab 



Ay^*.(r) 



ElV^,.(r) 



V(7(r)^(r)/5.(r; 



Ap(r) + Ap^(r; 



Pa(r)7(r)^(r) 



(26) 



where a denotes the spin projection opposite to a, i.e. 5" =1 if a =1 and vice 
versa. 

In the numerical calculations, equation (|T6|) is solved self-consistently together 
with the local one-particle equation 



1. 



: V + VI 



OEP( 



where 



V' 



OEP/ 



xccr V / J, 



r — r 



dr'. 



(27) 



(2^ 



3 Results 

Equation (pTf ) is solved by expanding the single-particle orbitals in spherical har- 
monics. The resulting equation for the radial part of the wave functions is then 
solved numerically. The main features of the code are as described in Ref. 28. 

For comparison, we have also performed calculations using the conventional 
Kohn-Sham method with two of the best standard exchange-correlation energy 
functionals. The first one of these is the exchange-energy functional by Becke 

combined with the correlation-energy functional by Lee, Yang and Parr 
referred to as BLYP. The other is the generalized gradient approximation by 
Perdew and Wang 0, referred to as PW91. 

Table |I] shows the total absolute ground-state energies of the first-row atoms. 
For these atoms, there exist accurate estimates of the exact non-relativistic values 
obtained from experimental ionisation energies and improved ab initio calcula- 
tions by Davidson et al. It is evident from the Table, that the density 



functional methods perform quite well. The mean absolute errors. A, given in 
the last row of Table |l|, clearly show that the present approach is significantly 
more accurate than the conventional Kohn-Sham method and nearly as accu- 



rate as recent CI results by Montgomery et al. The situation is similar 

for second-row atoms, as can be seen from Table 0. As the relativistic effects 
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Table 1: Total absolute ground-state energies for first-row atoms from various 
self-consistent calculations. CI values from Ref 29. A denotes the mean absolute 
deviation from the exact nonrelativistic values [30]. All numbers in Hartree units. 





present 


BLYP 


PW91 


CI 


EXACT 


He 


2.9033 


2.9071 


2.9000 


2.9049 


2.9037 


Li 


7.4829 


7.4827 


7.4742 


7.4743 


7.4781 


Be 


14.6651 


14.6615 


14.6479 


14.6657 


14.6674 


B 


24.6564 


24.6458 


24.6299 


24.6515 


24.6539 


C 


37.8490 


37.8430 


37.8265 


37.8421 


37.8450 


N 


54.5905 


54.5932 


54.5787 


54.5854 


54.5893 





75.0717 


75.0786 


75.0543 


75.0613 


75.067 


F 


99.7302 


99.7581 


99.7316 


99.7268 


99.734 


Nc 


128.9202 


128.9730 


128.9466 


128.9277 


128.939 


A 


0.0047 


0.0108 


0.0114 


0.0045 





Table 2: Total absolute ground-state energies for second-row atoms from various 
self-consistent calculations. A denotes the mean absolute deviation from Lamb- 
shift corrected experimental values, taken from Ref. 3. All numbers in Hartree 
units. 





present 


BLYP 


PW91 


EXPT 


Na 


162.256 


162.293 


162.265 


162.257 


Mg 


200.062 


200.093 


200.060 


200.059 


Al 


242.362 


242.380 


242.350 


242.356 


Si 


289.375 


289.388 


289.363 


289.374 


P 


341.272 


341.278 


341.261 


341.272 


S 


398.128 


398.128 


398.107 


398.139 


CI 


460.164 


460.165 


460.147 


460.196 


Ar 


527.553 


527.551 


527.539 


527.604 


A 


0.013 


0.026 


0.023 
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Table 3: Exchange and correlation energies from various approximations. All 
values in Hartree units. 





j^pvesent 










Tpt^w yi 


He 


1.028 


1.018 


1.009 


0.0416 


0.0437 


0.0450 


Li 


1.784 


1.771 


1.758 


0.0509 


0.0541 


0.0571 


Be 


2.674 


2.658 


2.644 


0.0934 


0.0954 


0.0942 


B 


3.760 


3.727 


3.711 


0.1289 


0.1287 


0.1270 


C 


5.064 


5.028 


5.010 


0.1608 


0.1614 


0.1614 


N 


6.610 


6.578 


6.558 


0.1879 


0.1925 


0.1968 





8.200 


8.154 


8.136 


0.2605 


0.2640 


0.2587 


F 


10.025 


9.989 


9.972 


0.3218 


0.3256 


0.3193 


Ne 


12.110 


12.099 


12.082 


0.3757 


0.3831 


0.3784 


Na 


14.017 


14.006 


13.985 


0.4005 


0.4097 


0.4040 


Mg 


15.997 


15.986 


15.967 


0.4523 


0.4611 


0.4486 


Al 


18.081 


18.053 


18.033 


0.4905 


0.4979 


0.4891 


Si 


20.295 


20.260 


20.238 


0.5265 


0.5334 


0.5322 


P 


22.649 


22.609 


22.587 


0.5594 


0.5676 


0.5762 


S 


25.021 


24.967 


24.944 


0.6287 


0.6358 


0.6413 


CI 


27.530 


27.476 


27.453 


0.6890 


0.6955 


0.7055 


Ar 


30.192 


30.139 


30.116 


0.7435 


0.7515 


0.7687 



for these atoms are more important and experiments increasingly difficult, the 
comparison of the calculated values with the Lamb-shift corrected experimental 
ones (from Ref. 3) has to be done cautiously and is by no means as rigorous 
as for ffist-row atoms. Nevertheless, our calculated values seem to mirror these 
experimental values more closely than the other approximations. 

For further analysis, we list, in Table ^, the values of and E^. separately. 
The data show two main features: First, the results for E^ are lowest for our 
method and highest for PW91, while the BLYP- values lie somewhere in between. 
And second, for E^., this trend is reversed, as now our results are highest and the 
ones from BLYP and PW91 are lower in nearly all cases. In Table ^ we show 
results of various x-only calculations performed with only the exchange-energy 
parts of the respective functionals. For the spherical atoms listed, there exist 



exact x-only OEP values [^, It is evident, that the KLI-approximation gives 
values much closer to the exact ones than the generalized gradient approxima- 
tions. From this and from Table ^ one may conclude that in the BLYP and PW91 
schemes an error cancellation between exchange and correlation energies occurs 
which leads to rather good total energies. Exchange and correlation energies 
separately, however, are poorly reproduced in the BLYP and PW91 approaches. 
By contrast, in our scheme, both exchange and correlation energies are of high 
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Table 4: Total absolute exchange energies of spherical first and second row atoms 
for various self consistent x-only calculations. The exact OEP data are from Refs. 
22 and 23. All values in Hartree units. 





KLI 


B88 


PW91 


OEP 


He 


1.026 


1.016 


1.005 


1.026 


Li 


1.781 


1.768 


1.754 


1.781 


Be 


2.667 


2.652 


2.638 


2.666 


N 


6.603 


6.569 


6.547 


6.604 


Ne 


12.099 


12.086 


12.061 


12.105 


Na 


14.006 


13.993 


13.968 


14.013 


Mg 


15.983 


15.972 


15.950 


15.988 


P 


22.633 


22.593 


22.565 


22.634 


Ar 


30.174 


30.122 


30.089 


30.175 



quality. 

Limitations of the three DFT approaches studied here become evident for 
ionisation potentials and electron affinities. In Table |^ we show ionisation poten- 
tials calculated from ground-state-energy differences and CI values from Ref. 29 
as well as experimental ones from Ref. 31. The performance of the three DFT 
methods is about the same, while the CI approach leads to clearly better results. 
Somewhat surprisingly, the DFT methods work better for the second-row than 
for the first-row atoms. 

In exact DFT, the highest occupied orbital energy of the neutral atom is iden- 
tical with the ionisation potential, while for negative ions the highest occupied 
energy level coincides with the electron affinity of the neutral atom |^ . How well 
ionisation potentials and electron affinities are reproduced by the highest occu- 
pied energy eigenvalues resulting from an approximate xc functional is therefore 
a measure of the quality of the functional. Table ^ shows the ionisation poten- 
tials obtained from the highest occupied single-particle-energy eigenvalue of the 
neutral atoms. The resulting values are worse than the ones in Table ^. While 
the deviation from experiment resulting from the BLYP and PW91 functionals 
is around 100 percent for all atoms, it is at most 10 percent in our approach. 
This is due to the fact that in our approach the xc potential has the correct 
— - asymptotic behaviour for large r which is not properly reproduced with the 
BLYP and PW91 functionals. 

For electron affinities, the situation is much worse, as may be seen from Table 
1^. First of all, because of the wrong asymptotic behaviour of the xc potential 
for large r, there is no convergence for negative ions within the self-consistent 
BLYP and PW91 schemes. This is not the case for our approach. However, the 
resulting electron affinities obtained either from ground-state-energy differences 
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Table 5: lonisation potentials calculated from ground-state-energy differences 
of neutral atoms. CI values are from Ref 29. A denotes the mean absolute 
deviation from the experimental values, taken from Ref. 31. All values in Hartree 
units. 





present 


BLYP 


PW91 


CI 


EXPT 


He 


0.903 


0.912 




0.905 


0.903 


Li 


0.203 


0.203 


0.207 


0.198 


0.198 


Be 


0.330 


0.330 


0.333 


0.344 


0.343 


B 


0.314 


0.309 


0.314 


0.304 


0.305 


C 


0.414 


0.425 


0.432 


0.413 


0.414 


N 


0.527 


0.542 


0.551 


0.534 


0.534 





0.495 


0.508 


0.505 


0.499 


0.500 


F 


0.621 


0.656 


0.660 


0.639 


0.640 


Ne 


0.767 


0.808 


0.812 


0.792 


0.792 


A 


0.009 


0.010 


0.014 


0.001 




Na 


0.191 


0.197 


0.198 




0.189 


Mg 


0.275 


0.280 


0.281 




0.281 


Al 


0.218 


0.212 


0.221 




0.220 


Si 


0.294 


0.294 


0.305 




0.300 


P 


0.379 


0.376 


0.389 




0.385 


S 


0.380 


0.379 


0.379 




0.381 


CI 


0.471 


0.476 


0.482 




0.477 


Ar 


0.575 


0.576 


0.583 




0.579 


A 


0.004 


0.005 


0.004 







or from the highest orbital energies of negative ions are not very reliable, in fact, 
for the Boron atom, the wrong sign is obtained from the ground-state-energy 
differences. Here, quantum-chemical approaches, such as CI, are clearly superior. 

4 Conclusions 

We have studied a novel density-functional scheme combining the KLI approxi- 
mation for the optimized effective potential with the correlation-energy functional 
of CoUe and Salvetti. The total ground-state energies obtained with this method 
for first and second-row atoms are significantly more accurate than the ones from 
the standard BLYP and PW91 schemes and only slightly less accurate than recent 
CI calculations. Error cancellations between exchange and correlation energies, 
well known to be present in the standard LDA, BLYP and PW91 functionals are 
not found for the new scheme. Therefore, both exchange and correlation energies 



11 



Tabic 6: lonisation potentials from the highest occupied orbital energy of neutral 
atoms. Experimental values from Ref. 31. All values in Hartree units. 





present 


BLYP 


PW91 


EXPT 


He 


0.945 


0.585 


0.583 


0.903 


Li 


0.200 


0.111 


0.119 


0.198 


Be 


0.329 


0.201 


0.207 


0.343 


B 


0.328 


0.143 


0.149 


0.305 


c 


0.448 


0.218 


0.226 


0.414 


N 


0.579 


0.297 


0.308 


0.534 





0.559 


0.266 


0.267 


0.500 


F 


0.714 


0.376 


0.379 


0.640 


Ne 


0.884 


0.491 


0.494 


0.792 


Na 


0.189 


0.106 


0.113 


0.189 


Mg 


0.273 


0.168 


0.174 


0.281 


Al 


0.222 


0.102 


0.112 


0.220 


Si 


0.306 


0.160 


0.171 


0.300 


P 


0.399 


0.219 


0.233 


0.385 


S 


0.404 


0.219 


0.222 


0.381 


CI 


0.506 


0.295 


0.301 


0.477 


Ar 


0.619 


0.373 


0.380 


0.579 



Table 7: Self-consistent electron affinities, (a) from ground-state-energy differ- 
ences, and (b ) from the highest occupied orbital energies of the negative ions. CI 
values are from Ref. 29, experimental ones from Ref. 31. All values in Hartree 
units. 





present" 


present 


CI 


EXPT 


Li 


0.016 


0.024 


0.023 


0.023 


B 


-0.002 


0.033 


0.008 


0.010 


C 


0.028 


0.083 


0.045 


0.046 





0.017 


0.110 


0.052 


0.054 


F 


0.082 


0.208 


0.125 


0.125 


Na 


0.015 


0.022 




0.020 


Al 


0.007 


0.024 




0.016 


Si 


0.040 


0.065 




0.051 


P 


0.022 


0.048 




0.027 


S 


0.065 


0.106 




0.076 


CI 


0.122 


0.174 




0.133 
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are separately of high accuracy. 

The calculated atomic ionisation potentials are in satisfactory agreement with 
experiment but clearly inferior to CI results. If the ionisation potentials are 
calculated by taking ground-state-energy differences, the performance of the three 
DFT schemes is about the same. If, on the other hand, the ionisation potentials 
are calculated from the highest occupied orbital energy the new scheme is clearly 
superior to the BLYP and PW91 schemes. Percentage deviations are better by 
roughly an order of magnitude. This can be explained by the fact that the KLI 
approximation preserves the correct asymptotic — - behaviour of the exact Kohn- 
Sham exchange potential while the BLYP and PW91 functionals are deficient in 
this respect. For the same reason, negative ions are not bound in the BLYP 
and PW91 self-consistent schemes so that electron affinities cannot be obtained, 
whereas the new scheme allows the calculation of these quantities. The calculated 
values, however, are not very reliable in comparison with experimental data and 
CI results. This clearly shows the need of further improving the correlation- 
energy functional. Work along these lines is in progress. 

We finally emphasize that the numerical effort involved in the proposed scheme 
is only marginally higher than that of an ordinary Hartree-Fock calculation, 
whereas all the principles of standard DFT are preserved. We therefore expect 
that our approach can be used with great success in chemical studies of more 
complex and larger systems, improving the results []TU|, |TI| obtained with the 
ordinary Kohn-Sham method. 



5 Acknowledgements 

We gratefully appreciate the help of Dr. E. Engel especially for providing us with 
a conventional Kohn-Sham computer code and for some helpful discussions. We 
would also like to thank Professor J. Perdew for providing us with the PW91 xc 
subroutine, and Professor H. StoU and M. Petersilka for many useful discussions. 
This work was supported in part by the Deutsche Forschungsgemeinschaft. 



References 

[1] p. Hohenberg, W. Kohn, Phys. Rev. 136, B864 (1964) 
[2] W. Kohn, L.J. Sham, Phys. Rev. 140, A1133 (1965) 

[3] R.M. Dreizler, E.K.U. Gross, Density Functional Theory: An Approaeh to 
the Quantum Many-Body Problem ( Springer- Verlag, Berlin, 1990) 

[4] R.G. Parr, W. Yang, Density- Functional Theory of Atoms and Molecules 
(Oxford University, New York, 1989) 



13 



[5] Density Functional Methods in Chemistry, edited by J.K. Labanowski and 
J.W. Andzelm ( Springer- Verlag, Berlin 1991) 

[6] S.J. Vosko, L. Wilk, M. Nusair, Can. J. Phys. 58, 1200 (1980) 

[7] J. P. Perdew in Electronic structure of solids '91. edited by P. Ziesche and H. 
Eschrig (Akademie Verlag, Berlin 1991); J. P. Predew and Y. Wang, unpub- 
lished. 

[8] A.D. Becke, Phys. Rev. A38, 3098 (1988) 
[9] C. Lee, W. Yang, R.G. Parr, Phys. Rev. B37, 785 (1988) 
[10] B.C. Johnson, P.M.W. Gill, J.A. Pople, J. Chem. Phys. 98, 5612 (1992) 
[11] P.M.W. Gill, B.G. Johnson, J.A. Pople, Chem. Phys. Lett. 197, 499 (1992) 
[12] R.T. Sharp, G.K. Horton, Phys. Rev. 90, 317, (1953) 
[13] J.D. Talman, W.F. Shadwick, Phys. Rev. A14, 36 (1976) 
[14] T.L. Gilbert, Phys. Rev. B12, 2111 (1975) 

[15] Density Matrices and Density Functionals, edited by R. Erdahl and V.H. 
Smith, Jr. (D. Reidel, Dordrecht 1987) 

[16] J.E. Harriman, Phys. Rev. A24, 680 (1981) 

[17] G. Zumbach, K. Maschke, Phys. Rev. A23, 544 (1983); Erratum: Phys. Rev. 
A29, 1585 (1984) 

[18] J.B. Krieger, Y. Li, G.J. lafrate, Phys. Rev. A45, 101 (1992) 

[19] J.B. Krieger, Y. Li, G.J. lafrate, Phys. Rev. A46, 5453 (1992) 

[20] J.B. Krieger, Y. Li, G.J. lafrate, Phys. Rev. A47, 165 (1993) 

[21] M.R. Norman, D.D. Koelling, Phys. Rev. B30, 5530 (1984) 

[22] E. Engel, J.A. Chevary, L.D. Macdonald, S.H. Vosko, Z. Phys. D 23, 7 (1992) 

[23] E. Engel, S.H. Vosko, Phys. Rev. A47, 2800 (1993) 

[24] R. CoUe, D. Salvetti, Theor. Chim. Acta 37, 329 (1975) 

[25] R. CoUe, D. Salvetti, Theor. Chim. Acta 53, 55 (1979) 

[26] T. Zhu, C. Lee, W. Yang, J. Chem. Phys. 98, 4814 (1993) 

[27] J.C. Slater, Phys. Rev. 81, 385 (1951) 

14 



[28] C. Proese Fischer, The Hariree-Fock method for Atoms (Wiley, New York, 
1977) 

[29] J. A. Montgomery, J.W. Ochterski, G.A. Petersson, J. Chem. Phys. 101, 
5900 (1994) 

[30] E.R. Davidson, S.A. Hagstrom, S.J. Chakravorty, V.M. Umar, C. Froese 
Fischer, Phys. Rev. A 44, 7071 (1991) 

[31] A. A. Radzig, B.M. Smirnov, Reference Data on Atoms and Molecules 
(Springer Verlag, Berhn, 1985) 

[32] CO. Almbladh, U. von Barth, Phys. Rev. B31, 3231 (1985) 



15 



